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COMPRESSION  BY  LINERS  OF  PLASMA-FIELD 
CONFIGURATIONS  WITH  (3  < 1 


I.  Introduction 

A simple  dynamic  model  of  a cylindrical,  incompressible  liner 
imploding  into  a P = 1 plasma  was  developed  early  in  the  LINUS  pro- 
gram1 in  order  to  derive  the  basic  scaling  laws  for  power-producing 
systems.  This  model  has  recently  been  further  developed  to  cover  com- 
pressible liners,2  but  the  assumption  of  a P = 1 plasma  was  retained. 
Since  the  thick,  captive  liquid  liners  which  now  form  the  basis  of  the 
LINUS  concept3  are  inevitably  rather  slow,  axial  confinement  of  the 
plasma  is  essential,  and  one  method  of  achieving  this  is  by  means  of 
a closed  field-line  system  such  as  the  reversed-field  theta-pinch  (or 
belt-pinch,  if  an  azimuthal  field  is  included).  In  this  system  mag- 
netic field  and  plasma  are  mixed,  that  is  to  say,  3 < 1;  the  most  im- 
portant consequence  is  that,  for  a given  Q,  the  radius  of  the  system 
will  be  greater  than  in  the  3=1  case.  Clearly  it  is  desirable  to 
choose  conditions  that  minimize  this  increase,  in  order  to  maintain 
our  claim  that  LINUS  reactors  will  be  very  compact  devices. 

The  principal  relationship  obtained  from  the  incompressible 
model  of  Ref.  1 was 

Q=Br  o"P(b,  T ) (1) 

o o o 

Note:  Manuscript  submitted  April  6,  1977. 
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where  Q = nuclear  energy  released  in  one  cycle/plasma  energy  at  peak 
compression; 

= buffer  magnetic  field  between  liner  and  plasma  at  peak  com- 
pression (assumed  of  negligible  extent); 

r = minimum  liner  radius; 
o 

p = density  of  liner  material; 

T = plasma  temperature  at  peak  compression; 
o 

b = a dimensionless  thickness  parameter ,' defined  by  A = rrb2r2  , 
where  A is  the  cross-sectional  area  of  the  liner. 

The  function  P(b,  T ) was  evaluated  by  integrating  the  thermo- 
nuclear reaction  rate  over  a compression-expansion  cycle;  the  plasma 
was  assumed  to  have  uniform  density  and  temperature. 

When  plasma  and  field  are  mixed  in  a 3 < 1 system,  the  basic  re- 
lationship (1)  still  holds,  but  the  value  of  P is  modified  as  the  result 
of  a number  of  effects,  namely: 

1)  The  average  plasma  density  is  reduced,  reducing  the  fusion 
reaction  rate  and  hence  decreasing  P. 

2)  The  system  energy,  for  a given  pressure,  is  reduced  since  now 
5/3  < Y < 2,  compared  to  Y = 5/3  for  the  3=1  case.  This 
increases  P. 

3)  A consequence  of  (2)  is  that  the  expansion  velocity  is 
reduced,  and  the  dwell  time  increased.  This  further  increases 
P. 

4)  The  increase  in  V causes  the  pressure  to  fall  more  rapidly  as 
the  system  expands.  This  decreases  P. 
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Of  these  effects,  (!)  is  dominant  and  results  in  a net  decrease  in 
P.  However  the  other  effects  are  all  significant  and  will  be  included 
in  the  following  discussion. 


II.  Compression  of  a Reversed-Field  Pinch 

The  reversed-field  theta-pinch  (Fig.  1)  is  a long  cylindrical  con- 
figuration of  radius  R in  which  the  magnetic  field  B is  in  the  z-direc- 
tion  over  most  of  the  length,  and  goes  from  -Bi  on  axis  to  +Bi  at 
r = R.  The  field  lines  close  at  the  ends,  and  for  the  purpose  of  this 
report  we  assume  that  this  provides  perfect  axial  confinement  of  plasma. 
The  radial  distribution  of  plasma  pressure  p = 2nkT  is  determined  by 
the  pressure  balance  equation 

p + B2/8tt  = Pi  = B2/8tt  (2) 

There  is  a cylindrical  neutral  surface  in  the  center  of  the  plasma  at 
which  B = 0 and  here  the  local  (3,  defined  by  (3  = 2nkT/pa  will  be  unity. 
Elsewhere  |3  < 1 , falling  to  zero  at  r = 0 and  r = R. 

Consider  the  system  initially  at  pressure  p^  and  take  an  element 
of  plasma  initially  having  f3  = and  volume  AV^.  Then  if  the  pressure 
is  increased  to  p the  volume  AV  will  be  given  by 


the  terms  on  the  right-hand  side  representing  the  partial  pressures  of 

the  plasma  and  field  respectively.  This  may  be  written  in  the  form: 

3/5  , , ,-1/3  ,3/5 


AV 

AV 


i -(-y  + ° ■ (%)  ) 


<M 
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which  allows  = f / , (3 . ] to  be  obtained  by  successive  approxi- 

i Vpi  V 

mat  ions . 

The  volume  of  the  whole  system  is  then  found  by  integration: 


To  determine  V as  a function  of  p requires  that  the  initial  distribu- 
tion of  0 is  known,  and  this  in  turn  depends  on  the  method  of  creating 
the  initial  plasma.  For  LINUS  systems,  the  rotating  relativistic 
electron  beam  seems  a promising  method  of  inducing  the  plasma  currents 
necessary  to  create  the  reversed-f ield  configuration, 3 but  at  present 
little  is  known  about  the  pressure  distribution  so  formed.  In  this 
report  we  shall  therefore  consider  two  simple  representative  distribu- 
tions, corresponding  respectively  to  a uniform  distribution  of  current 
density  and  a parabolic  distribution  centered  on  the  neutral  surface, 
which  is  taken  to  be  at  r = R/2.  The  two  distributions  are  then: 


(a*)  j.  = constant;  0.  - 1 - 4x2/R2 
3 1 

(bl  j 1 - 4x2/r2;  p = 1 - (5x/R  - 4x3/r3) 2 (6) 

W 1- 


where  x is  the  distance  from  the  neutral  surface. 

We  have  used  equations  (4)  and  (51  to  derive  V/V^  as  a function  of 
p/'p^  for  these  two  distributions.  One  significant  result  is  that  y,  the 
ratio  of  the  specific  heats  of  the  plasma-field  mixture,  varies  with  the 
degree  of  compression.  Here 


(lnp'> 

d(£nV) 
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which  is  plotted  against  V^/V,  the  volume  compression  ratio,  in  Fig.  2; 


this  result  will  be  used  later  in  determining  the  liner  trajectories. 

111.  Thermonuclear  Reactions 

We  wish  to  compare  the  total  fusion  reaction  rate  W integrated 
over  the  plasma  cross-section,  with  the  corresponding  rate  in  a p = 1 
plasma.  The  comparison  is  greatly  simplified  by  the  fact  that,  for 
the  D-T  reaction,  gv  oc  1"  , witli  an  accuracy  of  ± 101’,  for  8 keV  < T < 
keV.  Provided  that  the  final  compressed  plasma  temperature  is  in 
this  range,  the  local  reaction  rate  is  proportional  to  n2!2  = p and 
we  need  only  consider  the  radial  distribution  of  pressure , without 
requiring  that  the  density  and  temperature  be  specified  independently. 
Thus , 
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Using  the  method  of  Section  II,  we  have  calculated  (p2)  as  a func- 
tion of  V^/V  for  the  two  initial  distributions  used  previously.  The 
results  are  given  in  Fig.  5-  It  is  clear  that  there  is  a very  signif- 
icant reduction  in  the  fusion  rate  at  larger  compressions,  especially 
for  distribution  (M . 

In  Ref.  1 the  liner  trajectories  were  calculated  assuming  that  the 
liner  was  compressing  a p ^ 1 plasma  with  Y = 5/3*  We  now  generalize 
this  calculation  to  arbitrary  v,  but  retain  the  assumption  of  p - 1. 

The  normalized  equation  of  motion  of  the  liner  is  now 

’ n ( 1 ( b / a ' ' 


da  b / 1 - a~ ' ~ Y 
dT  3 Vn(  1 1 b /a')' 


who  re 


a = 


The  subscripts  o 


r/r  , t — tv  /r  and  v = (t.'p  /(V-l)b  p)^. 
o on  O a>  O r 

refer  to  values  at  peak  compression.  The  energy  per  unit  length  is 
given  by 

Eq  = nr,7p0/(V  - 1)  (10) 


anti  the  ratio  Q of  reaction  energy  released  in  one  cycle  to  is 
given  by 


B r o‘ 
o o 


- P(b,  TqI  V)  = 


vv-n 

32rr* 


T- 


i ^ dT 

a" 


cm 


where  E^  is  the  fusion  energy  per  reaction,  which  we  take  as  17.6  MeV. 
As  can  be  seen  from  equation  ('))  , the  integral  is  a function  of  v. 

Since  it  has  been  shown  previously1  that  P is  rather  insensitive  to  T 
(a  consequence  of  the  variation  of  ~ with  T mentioned  earlier)  we  take 

T = 20  keV  and  derive  P(V)  for  b — b , 8 and  11.  The  results  are 

o 

shown  in  Table  1.  We  find  that  the  ratio  sj  = P(Y)/P(3/3)  is  insensi- 
tive to  b,  and  so  it  is  plotted  as  a function  of  V only  in  Fig.  k. 

We  may  now  combine  the  results  of  Figs.  2,  8 anil  1*  to  obtain 


s 


r - st  ov"> 
6-1 


(12) 


that  is,  the  ratio  of  P for  the  reversed  field  pinch  to  its  correspond- 
ing value  for  a fl  - 1 system.  s is  plotted  as  a function  of  V^/V  in 
Fig. 

Discussion 

The  foregoing  analysis  lias  been  arranged  in  such  a way  as  to 
identify  the  separate  contributions  of  the  various  effects  listed  in 
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Section  I.  The  dominant  effect  is  clearly  the  (P2)  term  (Fig.  3) . The 
effect  of  increased  y appears  at  first  in  Equation  (11)  as  a (y  - l)3'2 
term  which,  in  the  most  favorable  case  (y  = 1.95,  Fig.  2)  leads  to  an 
increase  in  P of  1.7,  but  when  the  integral  is  included  this  is  reduced 
to  1.22  (Fig.  4).  The  outcome  of  this  exercise,  culminating  in  Fig.  5, 
is  to  allow  us  to  draw  three  general  cone lus ions r 

1)  Uniform  initial  current  distributions  are  to  be  preferred 
to  peaked  distributions. 

2)  Volumetric  compression  ratios  should  probably  not  exceed 
about  100-200. 

3)  For  scaling  studies  of  reversed-field  LINUS  systems,  it  is 
appropriate  to  take  s = 0.2  as  the  value  that  could  probably 
be  achieved  in  practice.  For  the  most  optimistic  projections, 
one  might  take  s = 0.25. 

The  most  important  consequence  of  (3  < 1 is  an  increase  in  the 
radial  dimensions  of  a LINUS  reactor.  For  fixed  Q,  the  increase  is 
by  a factor  s_1;  however,  as  shown  in  Ref.  2,  as  the  radius  is  in- 
creased, the  Q necessary  for  a self-sustaining  reactor  cycle  decreases, 
and  the  net  increase  in  size  is  only  by  a factor  s typically 

about  3- 
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Fig.  1 — Reversed-field  theta-pinch  configuration 
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